Gene co-expression patterns in Atlantic salmon adipose tissue provide a molecular link among seasonal changes, energy balance and age at maturity.

Sexual maturation in many fishes requires a major physiological change that involves a rapid transition between energy storage and usage. In Atlantic salmon, this transition for the initiation of maturation is tightly controlled by seasonality and requires a high-energy status. Lipid metabolism is at the heart of this transition since lipids are the main energy storing molecules. The balance between lipogenesis (lipid accumulation) and lipolysis (lipid use) determines energy status transitions. A genomic region containing a transcription co-factor of the Hippo pathway, vgll3, is the main determinant of maturation timing in Atlantic salmon. Interestingly, vgll3 acts as an inhibitor of adipogenesis in mice and its genotypes are potentially associated with seasonal heterochrony in lipid storage and usage in juvenile Atlantic salmon. Here, we explored changes in expression of more than 300 genes directly involved in the processes of adipogenesis, lipogenesis and lipolysis, as well as the Hippo pathway in the adipose tissue of immature and mature Atlantic salmon with distinct vgll3 genotypes. We found molecular evidence consistent with a scenario in which immature males with different vgll3 genotypes exhibit contrasting seasonal dynamics in their lipid profiles. We also identified components of the Hippo signalling pathway as potential major drivers of vgll3 genotype-specific differences in adipose tissue gene expression. This study demonstrates the importance of adipose gene expression patterns for directly linking environmental changes with energy balance and age at maturity through genetic factors bridging lipid metabolism, seasonality and sexual maturation.

Gene expression patterns associated with caudal fin shape in the cichlid Lamprologus tigripictilis.

Variation in fin shape is one of the most prominent features of morphological diversity among fish. Regulation of fin growth has mainly been studied in zebrafish, and it is not clear whether the molecular mechanisms underlying shape variation are equally diverse or rather conserved across species. In the present study, expression levels of 37 candidate genes were tested for association with fin shape in the cichlid fish Lamprologus tigripictilis. The tested genes included members of a fin shape-associated gene regulatory network identified in a previous study and novel candidates selected within this study. Using both intact and regenerating fin tissue, we tested for expression differences between the elongated and the short regions of the spade-shaped caudal fin and identified 20 genes and transcription factors (including angptl5, cd63, csrp1a, cx43, esco2, gbf1, and rbpj), whose expression patterns were consistent with a role in fin growth. Collated with available gene expression data of two other cichlid species, our study not only highlights several genes that were correlated with fin growth in all three species (e.g., angptl5, cd63, cx43, and mmp9), but also reveals species-specific gene expression and correlation patterns, which indicate considerable divergence in the regulatory mechanisms of fin growth across cichlids. Supplementary Information: The online version contains supplementary material available at 10.1007/s10750-022-05068-4.

Conserved Molecular Players Involved in Human Nose Morphogenesis Underlie Evolution of the Exaggerated Snout Phenotype in Cichlids.

Instances of repeated evolution of novel phenotypes can shed light on the conserved molecular mechanisms underlying morphological diversity. A rare example of an exaggerated soft tissue phenotype is the formation of a snout flap in fishes. This tissue flap develops from the upper lip and has evolved in one cichlid genus from Lake Malawi and one genus from Lake Tanganyika. To investigate the molecular basis of snout flap convergence, we used mRNA sequencing to compare two species with snout flap to their close relatives without snout flaps from each lake. Our analysis identified 201 genes that were repeatedly differentially expressed between species with and without snout flap in both lakes, suggesting shared pathways, even though the flaps serve different functions. Shared expressed genes are involved in proline and hydroxyproline metabolism, which have been linked to human skin and facial deformities. Additionally, we found enrichment for transcription factor binding sites at upstream regulatory sequences of differentially expressed genes. Among the enriched transcription factors were members of the FOX transcription factor family, especially foxf1 and foxa2, which showed an increased expression in the flapped snout. Both of these factors are linked to nose morphogenesis in mammals. We also found ap4 (tfap4), a transcription factor showing reduced expression in the flapped snout with an unknown role in craniofacial soft tissue development. As genes involved in cichlid snout flap development are associated with human midline facial dysmorphologies, our findings hint at the conservation of genes involved in midline patterning across distant evolutionary lineages of vertebrates, although further functional studies are required to confirm this.

A pituitary gene network linking vgll3 to regulators of sexual maturation in male Atlantic salmon.

Age at maturity is a key life history trait and a significant contributor to life history strategy variation. The maturation process is influenced by genetic and environmental factors, but specific causes of variation in maturation timing remain elusive. In many species, the increase in the regulatory gonadotropin-releasing hormone 1 (GnRH1) marks the onset of puberty. Atlantic salmon, however, lacks the gnrh1 gene, suggesting gnrh3 and/or other regulatory factors are involved in the maturation process. Earlier research in Atlantic salmon has found a strong association between alternative alleles of vgll3 and maturation timing. Recently we reported strong induction of gonadotropin genes (fshb and lhb) in the pituitary of Atlantic salmon homozygous for the early maturation allele (E) of vgll3. The induction of gonadotropins was accompanied by increased expression of their direct upstream regulators, c-jun and sf1 (nr5a1b) but the regulatory connection between vgll3 and these regulators has never been investigated in any organism. In this study, we investigated the potential regulatory connection between vgll3 genotypes and these regulators through a stepwise approach of identifying a gene regulatory network (GRN) containing c-jun and sf1, and transcription factor motif enrichment analysis. We found a GRN containing c-jun with predicted upstream regulators, e2f1, egr1, foxj1 and klf4, to be differentially expressed in the pituitary. Finally, we suggest a vgll3 and Hippo pathway -dependent model for transcriptional regulation of c-jun and sf1 in the pituitary, which may have broader implications across vertebrates.

Appetite regulating genes in zebrafish gut; a gene expression study.

The underlying molecular pathophysiology of feeding disorders, particularly in peripheral organs, is still largely unknown. A range of molecular factors encoded by appetite-regulating genes are already described to control feeding behaviour in the brain. However, the important role of the gastrointestinal tract in the regulation of appetite and feeding in connection to the brain has gained more attention in the recent years. An example of such inter-organ connection can be the signals mediated by leptin, a key regulator of body weight, food intake and metabolism, with conserved anorexigenic effects in vertebrates. Leptin signals functions through its receptor (lepr) in multiple organs, including the brain and the gastrointestinal tract. So far, the regulatory connections between leptin signal and other appetite-regulating genes remain unclear, particularly in the gastrointestinal system. In this study, we used a zebrafish mutant with impaired function of leptin receptor to explore gut expression patterns of appetite-regulating genes, under different feeding conditions (normal feeding, 7-day fasting, 2 and 6-hours refeeding). We provide evidence that most appetite-regulating genes are expressed in the zebrafish gut. On one hand, we did not observed significant differences in the expression of orexigenic genes (except for hcrt) after changes in the feeding condition. On the other hand, we found 8 anorexigenic genes in wild-types (cart2, cart3, dbi, oxt, nmu, nucb2a, pacap and pomc), as well as 4 genes in lepr mutants (cart3, kiss1, kiss1r and nucb2a), to be differentially expressed in the zebrafish gut after changes in feeding conditions. Most of these genes also showed significant differences in their expression between wild-type and lepr mutant. Finally, we observed that impaired leptin signalling influences potential regulatory connections between anorexigenic genes in zebrafish gut. Altogether, these transcriptional changes propose a potential role of leptin signal in the regulation of feeding through changes in expression of certain anorexigenic genes in the gastrointestinal tract of zebrafish.

Strong regulatory effects of vgll3 genotype on reproductive axis gene expression in juvenile male Atlantic salmon.

Age at maturity is a major contributor to the diversity of life history strategies in organisms. The process of maturation is influenced by both genetics and the environment, and includes changes in levels of sex hormones and behavior, but the specific factors leading to variation in maturation timing are poorly understood. gnrh1 regulates the transcription of gonadotropin genes at pubertal onset in many species, but this gene is lacking in certain teleost species including Atlantic salmon (Salmo salar), which raises the possibility of the involvement of other important regulatory factors during this process. Earlier research has reported a strong association of alternative alleles of the vgll3 gene with maturation timing in Atlantic salmon, suggesting it as a potential candidate regulating reproductive axis genes. Here, we investigated the expression of reproductive axis genes in one-year-old Atlantic salmon males with immature gonads and different vgll3 genotypes during the spawning period. We detected strong vgll3 genotype-dependent differential expression of reproductive axis genes (such as fshb, lhb, amh and igf3) tested in the pituitary, and testis. In addition, we observed differential expression of jun (ap1) and nr5a1b (sf1), potential upstream regulators of gonadotropins in the pituitary, as well as axin2, id3, insl3, itch, ptgs2a and ptger4b, the downstream targets of amh and igf3 in the testis. Hereby, we provide evidence of strong vgll3 genotype-dependent transcriptional regulation of reproductive axis genes prior to sexual maturation and suggest alternative models for distinct actions of vgll3 genotypes on the related molecular processes.

Expression variations in ectodysplasin-A gene (eda) may contribute to morphological divergence of scales in haplochromine cichlids.

BACKGROUND: Elasmoid scales are one of the most common dermal appendages and can be found in almost all species of bony fish differing greatly in their shape. Whilst the genetic underpinnings behind elasmoid scale development have been investigated, not much is known about the mechanisms involved in moulding of scales. To investigate the links between gene expression differences and morphological divergence, we inferred shape variation of scales from two different areas of the body (anterior and posterior) stemming from ten haplochromine cichlid species from different origins (Lake Tanganyika, Lake Malawi, Lake Victoria and riverine). Additionally, we investigated transcriptional differences of a set of genes known to be involved in scale development and morphogenesis in fish. RESULTS: We found that scales from the anterior and posterior part of the body strongly differ in their overall shape, and a separate look on scales from each body part revealed similar trajectories of shape differences considering the lake origin of single investigated species. Above all, nine as well as 11 out of 16 target genes showed expression differences between the lakes for the anterior and posterior dataset, respectively. Whereas in posterior scales four genes (dlx5, eda, rankl and shh) revealed significant correlations between expression and morphological differentiation, in anterior scales only one gene (eda) showed such a correlation. Furthermore, eda displayed the most significant expression difference between species of Lake Tanganyika and species of the other two younger lakes. Finally, we found genetic differences in downstream regions of eda gene (e.g., in the eda-tnfsf13b inter-genic region) that are associated with observed expression differences. This is reminiscent of a genetic difference in the eda-tnfsf13b inter-genic region which leads to gain or loss of armour plates in stickleback. CONCLUSION: These findings provide evidence for cross-species transcriptional differences of an important morphogenetic factor, eda, which is involved in formation of ectodermal appendages. These expression differences appeared to be associated with morphological differences observed in the scales of haplochromine cichlids indicating potential role of eda mediated signal in divergent scale morphogenesis in fish.

Impaired leptin signaling causes subfertility in female zebrafish.

Reproduction is an energetically costly event across vertebrates and tightly linked to nutritional status and energy reserves. In mammals, the hormone leptin is considered as a link between energy homeostasis and reproduction. However, its role in fish reproduction is still unclear. In this study, we investigated the possible role of leptin in the regulation of reproduction in zebrafish, using a loss of function leptin receptor (lepr) strain. Impaired leptin signaling resulted in severe reproductive deficiencies in female zebrafish. lepr mutant females laid significantly fewer eggs, with low fertilization rates compared to wild-type females. Folliculogenesis was not affected, but oocyte maturation and ovulation were disrupted in lepr mutants. Interestingly, the expression of luteinizing hormone beta (lhb) in the pituitary was significantly lower in mutant females. Analysis of candidate genes in the ovaries and isolated fully grown follicles revealed differential expression of genes involved in steroidogenesis, oocyte maturation and ovulation in the mutants, which are known to be regulated by LH signaling. Moreover, subfertility in lepr mutants could be partially restored by administration of human chorionic gonadotropin. In conclusion, our results show that leptin deficiency does not affect early stages of follicular development, but leptin might be essential in later steps, such as in oocyte maturation and ovulation. To our knowledge, this is the first time that leptin is associated to reproductive deficiencies in zebrafish.

The importance of alternative splicing in adaptive evolution.

Although alternative splicing is a ubiquitous co-transcriptional gene regulatory mechanism in plants, animals and fungi, its contribution to evolutionary transitions is understudied. Alternative splicing enables different mRNA isoforms to be generated from the same gene, expanding transcriptomic and thus proteomic diversity. While the role of gene expression variation in adaptive evolution is widely accepted, biologists still debate the functional impact of alternative isoforms on phenotype. In light of recent empirical research linking splice variation to ecological adaptations, we propose that alternative splicing is an important substrate for adaptive evolution and speciation, particularly at short timescales. In this article we synthesise what is known about the role of alternative splicing in adaptive evolution. We discuss the contribution of standing splice variation to phenotypic plasticity and how hybridisation can produce novel splice forms. Going forwards, we propose that alternative splicing be included as a standard analysis alongside gene expression analysis so we can better understand of how alternative splicing contributes to adaptive divergence at the micro- and macroevolutionary levels.

Transcriptional study reveals a potential leptin-dependent gene regulatory network in zebrafish brain.

The signal mediated by leptin hormone and its receptor is a major regulator of body weight, food intake and metabolism. In mammals and many teleost fish species, leptin has an anorexigenic role and inhibits food intake by influencing the appetite centres in the hypothalamus. However, the regulatory connections between leptin and downstream genes mediating its appetite-regulating effects are still not fully explored in teleost fish. In this study, we used a loss of function leptin receptor zebrafish mutant and real-time quantitative PCR to assess brain expression patterns of several previously identified anorexigenic genes downstream of leptin signal under different feeding conditions (normal feeding, 7-day fasting, 2 and 6-h refeeding). These downstream factors include members of cart genes, crhb and gnrh2, as well as selected genes co-expressed with them based on a zebrafish co-expression database. Here, we found a potential gene expression network (GRN) comprising the abovementioned genes by a stepwise approach of identifying co-expression modules and predicting their upstream regulators. Among the transcription factors (TFs) predicted as potential upstream regulators of this GRN, we found expression pattern of sp3a to be correlated with transcriptional changes of the downstream gene network. Interestingly, the expression and transcriptional activity of Sp3 orthologous gene in mammals have already been implicated to be under the influence of leptin signal. These findings suggest a potentially conserved regulatory connection between leptin and sp3a, which is predicted to act as a transcriptional driver of a downstream gene network in the zebrafish brain.

Transcriptomics unravels molecular players shaping dorsal lip hypertrophy in the vacuum cleaner cichlid, Gnathochromis permaxillaris.

BACKGROUND: Teleosts display a spectacular diversity of craniofacial adaptations that often mediates ecological specializations. A considerable amount of research has revealed molecular players underlying skeletal craniofacial morphologies, but less is known about soft craniofacial phenotypes. Here we focus on an example of lip hypertrophy in the benthivorous Lake Tangnayika cichlid, Gnathochromis permaxillaris, considered to be a morphological adaptation to extract invertebrates out of the uppermost layer of mud bottom. We investigate the molecular and regulatory basis of lip hypertrophy in G. permaxillaris using a comparative transcriptomic approach. RESULTS: We identified a gene regulatory network involved in tissue overgrowth and cellular hypertrophy, potentially associated with the formation of a locally restricted hypertrophic lip in a teleost fish species. Of particular interest were the increased expression level of apoda and fhl2, as well as reduced expression of cyp1a, gimap8, lama5 and rasal3, in the hypertrophic lip region which have been implicated in lip formation in other vertebrates. Among the predicted upstream transcription factors, we found reduced expression of foxp1 in the hypertrophic lip region, which is known to act as repressor of cell growth and proliferation, and its function has been associated with hypertrophy of upper lip in human. CONCLUSION: Our results provide a genetic foundation for future studies of molecular players shaping soft and exaggerated, but locally restricted, craniofacial morphological changes in fish and perhaps across vertebrates. In the future, we advocate integrating gene regulatory networks of various craniofacial phenotypes to understand how they collectively govern trophic and behavioural adaptations.

Gene co-expression network analysis reveals mechanisms underlying ozone-induced carbamazepine toxicity in zebrafish (Danio rerio) embryos.

Sewage effluent ozonation can reduce concentrations of chemical pollutants including pharmaceutical residues. However, the formation of potentially toxic ozonation byproducts (OBPs) is a matter of concern. This study sought to elucidate toxicity mechanisms of ozonated carbamazepine (CBZ), an anti-epileptic drug frequently detected in sewage effluents and surface water, in zebrafish embryos (Danio rerio). Embryos were exposed to ozonated and non-ozonated CBZ from 3 h post-fertilization (hpf) until 144 hpf. Embryotoxicity endpoints (proportion of dead and malformed embryos) were assessed at 24, 48, and 144 hpf. Heart rate was recorded at 48 hpf. Exposure to ozonated CBZ gave rise to cardiovascular-related malformations and reduced heart rate. Moreover, embryo-larvae exposed to ozonated CBZ displayed a lack of swim bladder inflation. Hence, the expression patterns of CBZ target genes involved in cardiovascular and embryonal development were investigated through a stepwise gene co-expression analysis approach. Two co-expression networks and their upstream transcription regulators were identified, offering mechanistic explanations for the observed toxicity phenotypes. The study presents a novel application of gene co-expression analysis elucidating potential toxicity mechanisms of an ozonated pharmaceutical with environmental relevance. The resulting data was used to establish a putative adverse outcome pathway (AOP).

Gene coexpression networks reveal molecular interactions underlying cichlid jaw modularity.

BACKGROUND: The oral and pharyngeal jaw of cichlid fishes are a classic example of evolutionary modularity as their functional decoupling boosted trophic diversification and contributed to the success of cichlid adaptive radiations. Most studies until now have focused on the functional, morphological, or genetic aspects of cichlid jaw modularity. Here we extend this concept to include transcriptional modularity by sequencing whole transcriptomes of the two jaws and comparing their gene coexpression networks. RESULTS: We show that transcriptional decoupling of gene expression underlies the functional decoupling of cichlid oral and pharyngeal jaw apparatus and the two units are evolving independently in recently diverged cichlid species from Lake Tanganyika. Oral and pharyngeal jaw coexpression networks reflect the common origin of the jaw regulatory program as there is high preservation of gene coexpression modules between the two sets of jaws. However, there is substantial rewiring of genetic architecture within those modules. We define a global jaw coexpression network and highlight jaw-specific and species-specific modules within it. Furthermore, we annotate a comprehensive in silico gene regulatory network linking the Wnt and AHR signalling pathways to jaw morphogenesis and response to environmental cues, respectively. Components of these pathways are significantly differentially expressed between the oral and pharyngeal jaw apparatus. CONCLUSION: This study describes the concerted expression of many genes in cichlid oral and pharyngeal jaw apparatus at the onset of the independent life of cichlid fishes. Our findings suggest that - on the basis of an ancestral gill arch network-transcriptional rewiring may have driven the modular evolution of the oral and pharyngeal jaws, highlighting the evolutionary significance of gene network reuse. The gene coexpression and in silico regulatory networks presented here are intended as resource for future studies on the genetics of vertebrate jaw morphogenesis and trophic adaptation.

Expression levels of the tetratricopeptide repeat protein gene ttc39b covary with carotenoid-based skin colour in cichlid fish.

Carotenoid pigments play a major role in animal body colouration, generating strong interest in the genes involved in the metabolic processes that lead from their dietary uptake to their storage in the integument. Here, we used RNA sequencing (RNA-Seq) to test for differentially expressed genes in a taxonomically replicated design using three pairs of related cichlid fish taxa from the genera Tropheus and Aulonocara. Within each pair, taxa differed in terms of red and yellow body colouration, and high-performance liquid chromatography (HPLC) analyses of skin extracts revealed different carotenoid profiles and concentrations across the studied taxa. Five genes were differentially expressed in all three yellow-red skin contrasts (dhrsx, nlrc3, tcaf2, urah and ttc39b), but only the tetratricopeptide repeat protein-coding gene ttc39b, whose gene product is linked to mammalian lipid metabolism, was consistently expressed more highly in the red skin samples. The RNA-Seq results were confirmed by quantitative PCR. We propose ttc39b as a compelling candidate gene for variation in animal carotenoid colouration. Since differential expression of ttc39b was correlated with the presence/absence of yellow carotenoids in a previous study, we suggest that ttc39b is more likely associated with the concentration of total carotenoids than with the metabolic formation of red carotenoids.

Comparative transcriptomics reveals candidate carotenoid color genes in an East African cichlid fish.

BACKGROUND: Carotenoids contribute significantly to animal body coloration, including the spectacular color pattern diversity among fishes. Fish, as other animals, derive carotenoids from their diet. Following uptake, transport and metabolic conversion, carotenoids allocated to body coloration are deposited in the chromatophore cells of the integument. The genes involved in these processes are largely unknown. Using RNA-Sequencing, we tested for differential gene expression between carotenoid-colored and white skin regions of a cichlid fish, Tropheus duboisi "Maswa", to identify genes associated with carotenoid-based integumentary coloration. To control for positional gene expression differences that were independent of the presence/absence of carotenoid coloration, we conducted the same analyses in a closely related population, in which both body regions are white. RESULTS: A larger number of genes (n = 50) showed higher expression in the yellow compared to the white skin tissue than vice versa (n = 9). Of particular interest was the elevated expression level of bco2a in the white skin samples, as the enzyme encoded by this gene catalyzes the cleavage of carotenoids into colorless derivatives. The set of genes with higher expression levels in the yellow region included genes involved in xanthophore formation (e.g., pax7 and sox10), intracellular pigment mobilization (e.g., tubb, vim, kif5b), as well as uptake (e.g., scarb1) and storage (e.g., plin6) of carotenoids, and metabolic conversion of lipids and retinoids (e.g., dgat2, pnpla2, akr1b1, dhrs). Triglyceride concentrations were similar in the yellow and white skin regions. Extracts of integumentary carotenoids contained zeaxanthin, lutein and beta-cryptoxanthin as well as unidentified carotenoid structures. CONCLUSION: Our results suggest a role of carotenoid cleavage by Bco2 in fish integumentary coloration, analogous to previous findings in birds. The elevated expression of genes in carotenoid-rich skin regions with functions in retinol and lipid metabolism supports hypotheses concerning analogies and shared mechanisms between these metabolic pathways. Overlaps in the sets of differentially expressed genes (including dgat2, bscl2, faxdc2 and retsatl) between the present study and previous, comparable studies in other fish species provide useful hints to potential carotenoid color candidate genes.

Transcriptional study of appetite regulating genes in the brain of zebrafish (Danio rerio) with impaired leptin signalling.

The hormone leptin is a key regulator of body weight, food intake and metabolism. In mammals, leptin acts as an anorexigen and inhibits food intake centrally by affecting the appetite centres in the hypothalamus. In teleost fish, the regulatory connections between leptin and other appetite-regulating genes are largely unknown. In the present study, we used a zebrafish mutant with a loss of function leptin receptor to investigate brain expression patterns of 12 orexigenic and 24 anorexigenic genes under different feeding conditions (normal feeding, 7-day fasting, 2 and 6-hours refeeding). Expression patterns were compared to wild-type zebrafish, in order to identify leptin-dependent differentially expressed genes under different feeding conditions. We provide evidence that the transcription of certain orexigenic and anorexigenic genes is influenced by leptin signalling in the zebrafish brain. We found that the expression of orexigenic genes was not affected by impaired leptin signalling under normal feeding conditions; however, several orexigenic genes showed increased transcription during fasting and refeeding, including agrp, apln, galr1a and cnr1. This suggests an inhibitory effect of leptin signal on the transcription of these orexigenic genes during short-term fasting and refeeding in functional zebrafish. Most pronounced effects were observed in the group of anorexigenic genes, where the impairment of leptin signalling resulted in reduced gene expression in several genes, including cart family, crhb, gnrh2, mc4r, pomc and spx, in the control group. This suggests a stimulatory effect of leptin signal on the transcription of these anorexigenic genes under normal feeding condition. In addition, we found multiple gain and loss in expression correlations between the appetite-regulating genes, in zebrafish with impaired leptin signal, suggesting the presence of gene regulatory networks downstream of leptin signal in zebrafish brain. The results provide the first evidence for the effects of leptin signal on the transcription of various appetite-regulating genes in zebrafish brain, under different feeding conditions. Altogether, these transcriptional changes suggest an anorexigenic role for leptin signal, which is likely to be mediated through distinct set of appetite-regulating genes under different feeding conditions.

Divergence in larval jaw gene expression reflects differential trophic adaptation in haplochromine cichlids prior to foraging.

BACKGROUND: Understanding how variation in gene expression contributes to morphological diversity is a major goal in evolutionary biology. Cichlid fishes from the East African Great lakes exhibit striking diversity in trophic adaptations predicated on the functional modularity of their two sets of jaws (oral and pharyngeal). However, the transcriptional basis of this modularity is not so well understood, as no studies thus far have directly compared the expression of genes in the oral and pharyngeal jaws. Nor is it well understood how gene expression may have contributed to the parallel evolution of trophic morphologies across the replicate cichlid adaptive radiations in Lake Tanganyika, Malawi and Victoria. RESULTS: We set out to investigate the role of gene expression divergence in cichlid fishes from these three lakes adapted to herbivorous and carnivorous trophic niches. We focused on the development stage prior to the onset of exogenous feeding that is critical for understanding patterns of gene expression after oral and pharyngeal jaw skeletogenesis, anticipating environmental cues. This framework permitted us for the first time to test for signatures of gene expression underlying jaw modularity in convergent eco-morphologies across three independent adaptive radiations. We validated a set of reference genes, with stable expression between the two jaw types and across species, which can be important for future studies of gene expression in cichlid jaws. Next we found evidence of modular and non-modular gene expression between the two jaws, across different trophic niches and lakes. For instance, prdm1a, a skeletogenic gene with modular anterior-posterior expression, displayed higher pharyngeal jaw expression and modular expression pattern only in carnivorous species. Furthermore, we found the expression of genes in cichlids jaws from the youngest Lake Victoria to exhibit low modularity compared to the older lakes. CONCLUSION: Overall, our results provide cross-species transcriptional comparisons of modularly-regulated skeletogenic genes in the two jaw types, implicating expression differences which might contribute to the formation of divergent trophic morphologies at the stage of larval independence prior to foraging.

Gene expression profiling suggests differences in molecular mechanisms of fin elongation between cichlid species.

Comparative analyses of gene regulation inform about the molecular basis of phenotypic trait evolution. Here, we address a fin shape phenotype that evolved multiple times independently across teleost fish, including several species within the family Cichlidae. In a previous study, we proposed a gene regulatory network (GRN) involved in the formation and regeneration of conspicuous filamentous elongations adorning the unpaired fins of the Neolamprologus brichardi. Here, we tested the members of this network in the blockhead cichlid, Steatocranus casuarius, which displays conspicuously elongated dorsal and moderately elongated anal fins. Our study provided evidence for differences in the anatomy of fin elongation and suggested gene regulatory divergence between the two cichlid species. Only a subset of the 20 genes tested in S. casuarius showed the qPCR expression patterns predicted from the GRN identified in N. brichardi, and several of the gene-by-gene expression correlations differed between the two cichlid species. In comparison to N. brichardi, gene expression patterns in S. casuarius were in better (but not full) agreement with gene regulatory interactions inferred in zebrafish. Within S. casuarius, the dorsoventral asymmetry in ornament expression was accompanied by differences in gene expression patterns, including potential regulatory differentiation, between the anal and dorsal fin.

Molecular mechanisms underlying nuchal hump formation in dolphin cichlid, Cyrtocara moorii.

East African cichlid fishes represent a model to tackle adaptive changes and their connection to rapid speciation and ecological distinction. In comparison to bony craniofacial tissues, adaptive morphogenesis of soft tissues has been rarely addressed, particularly at the molecular level. The nuchal hump in cichlids fishes is one such soft-tissue and exaggerated trait that is hypothesized to play an innovative role in the adaptive radiation of cichlids fishes. It has also evolved in parallel across lakes in East Africa and Central America. Using gene expression profiling, we identified and validated a set of genes involved in nuchal hump formation in the Lake Malawi dolphin cichlid, Cyrtocara moorii. In particular, we found genes differentially expressed in the nuchal hump, which are involved in controlling cell proliferation (btg3, fosl1a and pdgfrb), cell growth (dlk1), craniofacial morphogenesis (dlx5a, mycn and tcf12), as well as regulators of growth-related signals (dpt, pappa and socs2). This is the first study to identify the set of genes associated with nuchal hump formation in cichlids. Given that the hump is a trait that evolved repeatedly in several African and American cichlid lineages, it would be interesting to see if the molecular pathways and genes triggering hump formation follow a common genetic track or if the trait evolved in parallel, with distinct mechanisms, in other cichlid adaptive radiations and even in other teleost fishes.

Maternal mRNA input of growth and stress-response-related genes in cichlids in relation to egg size and trophic specialization.

BACKGROUND: Egg size represents an important form of maternal effect determined by a complex interplay of long-term adaptation and short-term plasticity balancing egg size with brood size. Haplochromine cichlids are maternal mouthbrooders showing differential parental investment in different species, manifested in great variation in egg size, brood size and duration of maternal care. Little is known about maternally determined molecular characters of eggs in fishes and their relation to egg size and trophic specialization. Here we investigate maternal mRNA inputs of selected growth- and stress-related genes in eggs of mouthbrooding cichlid fishes adapted to different trophic niches from Lake Tanganyika, Lake Malawi, Lake Victoria and compare them to their riverine allies. RESULTS: We first identified two reference genes, atf7ip and mid1ip1, to be suitable for cross-species quantification of mRNA abundance via qRT-PCR in the cichlid eggs. Using these reference genes, we found substantial variation in maternal mRNA input for a set of candidate genes related to growth and stress response across species and lakes. We observed negative correlation of mRNA abundance between two of growth hormone receptor paralogs (ghr1 and ghr2) across all haplochromine cichlid species which also differentiate the species in the two younger lakes, Malawi and Lake Victoria, from those in Lake Tanganyika and ancestral riverine species. Furthermore, we found correlations between egg size and maternal mRNA abundance of two growth-related genes igf2 and ghr2 across the haplochromine cichlids as well as distinct clustering of the species based on their trophic specialization using maternal mRNA abundance of five genes (ghr1, ghr2, igf2, gr and sgk1). CONCLUSIONS: These findings indicate that variations in egg size in closely related cichlid species can be linked to differences in maternal RNA deposition of key growth-related genes. In addition, the cichlid species with contrasting trophic specialization deposit different levels of maternal mRNAs in their eggs for particular growth-related genes; however, it is unclear whether such differences contribute to differential morphogenesis at later stages of development. Our results provide first insights into this aspect of gene activation, as a basis for future studies targeting their role during ecomorphological specialization and adaptive radiation.